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ABSTRACT 

Electrostatic emissions are often observed by spacecraft in 
the outer magnetosphere at frequencies between low electron gyro- 
harmonies. The emissions include the well-lmown odd half -harmonic 
emissions as well as weaker emissions which occur near the electron 
plasma frequency. We have constructed a scheme for classifying the 
emissions and have shown that a theoretical model which has been 
previously used to explain the "3/2" emissions can be extended to the 
other classes of emissions. All of the emissions appear therefore 
to be generated by the same basic mechanism: an unstable electron 
plasma distribution consisting of cold electrons (<100 eV) and hot 
loss cone electrons (~1 keV). Each emission class is associated with 
a particular range of model parameters; the wide band electric field 
data can thus be used to infer the density and temperature of the cold 
plasma component. The model predicts that gyroharmonic emissions 
near the plasma frequency require large cold plasma densities. 
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ELECTROSTATIC EMISSIONS BETWEEN ELECTRON 


GYROHARMONICS IN TRE OUTER MAGNETOSPHERE 
L INTRODUCTION 

Banded electrostatic emissions at frequencies above the local electron 
gyrofrequency are a prominent and interesting feature of the outer magneto- 
sphere. By far the most common of these emissions are those occurring be- 
tween fg and 2fg , These were first observed by the plasma wave experiment on 
Ogo 5 (Kennel, ^ al. , 1970; Fredricks and Scarf, 1973) as narrowband 

~ 0.1) , large amplitude (1-10 inV/m) waves occurring over the approxi- 
mate range 1.25 fg < f< 1.75 f^. Because this frequency interval is centered 
at 1.5 fg, these waves are often loosely termed"3/2” emissions. Similar 
"3/2” emissions have also been detected by plasma wave experiments on 
(Anderson and Maeda, 1977) and Imp 6 (Shaw and Gurnett, 1975). In the Imp 6 
experiment "3/2" emissions are sporadically observed on nearly every traver- 
sal of the outer magnetosphere; they are indeed also a common magnetospheric 
feature. It has been suggested that these waves play a dominant role in the 
precipitation of 1-10 keV electrons (Lyons, 1974; Ashour-Abdalla and Kennel, 
1976). 

There are several other types of electrostatic turbulence whose charac- 
teristics suggest that they are related to "3/2" emissions. There are first of 
all "n + 1/2" or "odd half harmonic" emissions with n integral and > i. Again 
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the coinage is loose: these are narrow band emissions occurring between nfg 
and (n + but significantly away from the gyroharmonics themselves. The 
Ogo 5 experiment occasionally observed "5/2” and "7/2" emissions, and 
{n 4- l/2) emissions with n as large as 8 are common in the Imp 6 data, though 
far less so than the "3/2" emissions. 

When "n + l/2" emissions are seen in the magnetosphere, all half har- 
monics n = 1, 2, ... N are usually sequentially present. However, Shaw and 
Gurnett (1975) report narrow band emissions between higher gyroharmonics at 
times when the "3/2" and other low half harmonics are absent. In such cases, 
the highest frequency band is often near the local electron plasma frequency fp, 
the latter being determined from the observed lower cutoff of nonthermal con- 
tinuum radiation (Gurnett and Shaw, 1973). 

Somewhat different from the narrow banded emissions discussed so far 
are "diffuse" band emissions reported by Shaw and Gurnett (1975). These 
emissions are also confined to bands between the low electron gyroharmonics 
and tend to occur just outside the plasmapause. They may have substantially 
wider bandwidths (Af<0.5 fg), and much lower intensities than do narrow 
band emissions. The highest observed frequency of diffuse band emissions is 
frequently near f p. 

Banded electromagnetic emissions near fg/2 have also been reported 
(Burtis and Helliwell, 1969; Maeda, , 1976). These emissions, which are 

presumably whistler waves, are seldom observed unless "3/2" emissions are 
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siinultaneously present, and a causal link between the two may exist. We shall 
not be concerned, however, with electromagnetic emissions in this paper. 

Part of our task has been to develop a classification scheme for the emis- 
sions discussed in the previous paragraphs. Five separate categories of waves 
have been identified. More importantly, we have developed a parameterized 
model of a linearly unstable plasma. By suitable choices of the parameters 
in the model, it is possible to duplicate in a more or less unique fashion the 
emissions in each category. 

The plasma model is essentially the one investigated previously by several 
other authors (Fredricks, 1971; Young, et , 1973; Young, 1975; Ashour- 
Abdalla, _et 1975; Ashour-Abdalla and Kennel, 1976). Because of the high 
frequency of the waves involved, we consider only electron dynamics. There 
are two electron populations; a "hot” (~1 keV) loss cone component, with den- 
sity Hjj and parallel and perpendicular temperatures T|j^ and Tj^j^, is the source 
of free energy for the instability. In addition, there are "cold" {<100 eV) elec- 
trons of density and temperature Tj,. The cold electrons, which originate 
from the ionosphere, facilitate the propagation of the instability. As might be 
expected for electrostatic emissions whose frequency is related to f^ , the waves 
which we consider propagate nearly perpendicular to B , i.e. the perpendicular 
and parallel wavenumbers kjL and k|[ satisfy kj^ » k|]. 

In this paper, we extend the linear analysis to investigate ^t 5 q)es of 
narrow band electrostatic turbulence observed in the outer magnetosphere. We 
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feel that diffuse band emissions can also be explained by this model, but our 
work here is as yet tentative and will be reported later. The theoretical analy- 
sis of narrowband emissions is coupled to specific Imp 6 electric field observa- 
tions. The plasma conditions which give rise to emissions in each of the pro- 
posed categories are theoretically investigated. We summarize the effects of 
varying each of the many free parameters in the problem on the morphology of 
each category of emission. This leads us to a 'TDest estimate” as to local plasma 
conditions when a given t 5 ^e of emission occurs and hence a predictive or diag- 
nostic capability. In particular, we are able to estimate the local cold plasma 
density, a parameter which is difficult to measure by direct plasma observations. 
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II, TYPICAL ELECTRIC FIELD OBSERVATIONS 

The essential features of "half -harmonic" electrostatic emissions are 
most clearly seen in wide band spectrograms of the t 3 ^e shown in Figure 1. 

The figure shows time compressed electric field data from the University of 
Iowa plasma wave experiment on Imp 6 during an outbound orbital pass. Simul- 
taneous magnetic field measurements from the NASA/GSFC magnetometer allow 
us to calculate fg. The most persistent feature is the narrow bandwidth emis- 
sions between f^ and 2f^ which are almost continously observed from 0500 to 
0800. During this time, the relative frequency [i = f/f^ of this "3/2" emission 
varies between 1. 2 and 1.45 as the g 3 mofrequency decreases from 4. 6 kHz to 
0. 92 kHz. Figure 1, however, shows only the last part of this time interval. 

The "3/2" emissions in Figure 1 are almost constant in frequency and are the 
dark trace occurring between 1 kHz and 2 kHz. 

Higher odd half harmonic emissions occur sporadically after 0665 with 
as many as seven bands simultaneously present. The precise relative frequen- 
cies of the emissions at 0715 — f/fg=jti - 1.4, 2.3, 3. 2, 4.2, 5.2, and 6.1 — 
are indicative of the looseness of the "odd half-harmonic" terminology. The 
emissions below 1 kHz are the concurrent electromagnetic emissions near 

fj/2. 

Figure 2 shows one of the types of emissions reported by Shaw and Gurnett 
(1975) . The sporadic narrow band emissions which first appear about 0525 are 
the main feature of interest. The shaded region above ~6 kHz is the nonthermal 
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contiinamn radiation (Gurnett and Shaw, 1973). This radiation is electromag- 
netic and is density trapped between the plasmapause and the magnetopause. Its 
lower cutoff is the local electron plasma frequency fp. The narrow band emis- 
sions here are well above f„. The three bands at 0540 are just below 7f„ , 8f 
and 9fg. Although in this case, the frequencies of the narrow band emissions are 
just below fp , similar bands are sometimes observed above the lower continuum 
cutoff. 

Although the continuum radiation is usually observed whenever narrow 
band emissions of the type shown in Figure 2 are present, the continuum emis- 
sions are sometunes missing (e. g. , Shaw and Gurnett, 1975, Fig. 4). This may 
be due to having fp exceed the upper frequency limit (10 kHz or 30 kHz) of the 
wide band spectrometer, but it is more likely that the narrow band emissions 
are occurring near fp, and that the continuum is simplj^ too weak to be observed 
(R. Shaw, private communication) . 

The presence of continuum radiation above fp is not correlated with the 
occurrence of electrostatic emissions. Continuum radiation may accompany 
any tj^pe of narrow bandwidth electrostatic emissions, diffuse emissions (e.g. , 
Shaw and Gurnett, 1975, Figs. 6 and 11), or it may even occur in the complete 
absence of electrostatic turbulence. These pro"*erties, plus the fact that the 
continuum radiation is electromagnetic in nature, suggest that the mechanisms 
for generating continuum and odd half -harmonic emissions are unrelated. 
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Based on the work of Shaw and Gurnett (1975) and onr own survey of 
Imp 6 0-10 kHz broadband data, we propose the classification scheme shown 
in Table 1. With the exception of Classes 1 and 2, the distinction between 
separate categories is an observationally natural one, and as we shall see, each 
class corresponds to a different configuration of our plasma model. The origi- 
nal Ogo 5 observations were primarily low latitude measurements and found 
typical frequencies f 1. 5 f g, i. e. , Class 2 emissions. Our analyses of sev- 
eral orbits of Imp 6 data, however, suggest that <f <2f^ waves are concen- 
trated below 1. 2 fgt for example, 60% of the narrow band emissions observed 
in Orbits 178-183 were Class 1 as opposed to Class 2 emissions. Our theore- 
tical work suggests that this difference in frequency distributions between the 
Ogo 5 and Imp 6 observations is a latitudinal effect. At the same radial distance 
from the Earth, the ratio fp/fg> upon which the frequency f of narrow emissions 
depends strongly, is smaller at the higher latitudes traversed by Imp 6 (if f p 
is constant or decreases with latitude). Thus, the plasma characteristics of 
Class 1 and Class 2 events are also different. 

Class 3 emissions are multiple ”(n + l/2)" emissions of the type shown 
in Figure 1. Class 4 emissions are the diffuse electrostatic emissions described 
by Shaw and Gurnett (1975). Finally, Class 5 emissions are narrow band emis- 
sions such as those in Figure 2 between cyclotron harmonics near fp. 

We have assumed that all Class 5 events have frequencies controlled by 
cyclotron harmonics even when only a single band near f p is observed. It is 
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possible that some of these emissions are in fact electron plasma oscillations. 
We shall see that it is theoretically possible to produce g 3 a'oharmonic emissions 
which occur in a narrow band right at f^. In principle one should be able to dis- 
tinguish between electron plasma oscillations and gyroharmonic emissions from 
the location in time of intensity nulls associated with the spin period of the 
spacecraft since E is predominantly perpendicular to B for gyroharmonic emis-- 

“ 3 “ 

sions and parallel to B for electron plasma oscillations. Unfortunately, the 
Glass 5 emissions don't usually show spin modulation (Shaw and Gurnett, 1975). 
Emissions occurring in two or more bands near such as Figure 2, can 
safely be considered gyroharmonic emissions, while emissions which occur in 
a single narrow band near the continuum cutoff could be electron plasma oscilla- 
tions or gyroharmonic emissions. It may be possible to distinguish the two by 
analyzing electron particle spectra. 
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TABLE 1 


Glassification Scheme for Electrostatic Emissions in the Outer Magnetosphere 


Class 

Frequency Range 

Comments 

1. Low "3/2” 

Kf/fg<1.2, Af/fg<0.1 

Most common emission class for 20^ < |A | ;$ 55°. 

2. “3/2” 

1.2<f/fg<2. Af/fg<0.1 

Most widely studied emission class. Probably 
most common emission class near equator. 

3. “n + 1/2” or multiple 
half harmonic 

f /f >2 

max' g 

Af/fg < 0.1 in each band 

Narrow bandwidth emissions in each gyrohar- 
monic band above fg up to the band (e.g.. 

Fig. 1). 

4. Diffuse band 

f /f <4 

max^ g ~ 

Af/fg ~ 0.5 in lowest band 

Less sharply defined emissions in the lowest 
(1 < f/fg < 2) band. May have emissions in one 
or two other bands. If continuum observed, 
f ~ f . 

max p 

5. “f ~ fp ” narrow band 

f f ~ f 

g max p 

Af/fg SO.l 

One or more bands near f^ are enhanced. Lower 
(e.g., “3/2”) bands usually missing. 



IIL THEORY 


The dispersion relation for electrostatic waves in a magnetized electron 
plasma, where ions form only a fixed neutralizing background, is in cylindrical 
coordinates, with the z-axis along B (Harris, 1959), 


0= 1 -■ 
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Here co = C47rne^ and = 1 e|B/mc = 2nf„ are the angular electron plasma 
and gyrofrequencies respectively and the electron distribution f( Uj^ , U|j ) is 
normalized to unity 

CO 

0 

Two different classes of f (Vj^, V||) which lead to unstable solutions (7 = 

Im CO > 0) to Eq. (1) with kj^ » k|| have been widely studied. The first type is 
characterized by an anisotropy in which the temperature Tj^ associated with 
the electronic motion perpendicular to B exceeds the equivalent parallel tempera- 
ture T||. While > 2T|| is suffic lent for instability in the ”3/2” band, the 
instability is restricted to the region 3/2 f^ < f < 2fg (Young et al. , 1973), and 
in order to account for observations below 3/2 f a complicated non-linear 

o 

coupling must be invoked (Oya, 1972; 1975), 

Instability at higher odd half harmonics requires a still more pronounced 
temperature anisotropy. It is primarily because such anisotropic distributions 


/^"ii 


dUjjfCuj^, U||) = 1 
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are not observed in the magnetosphere that we analyze this type of distribution 
function no further. 

The second category of distribution functions is one in which 3f/3vj^ > 0 
over some region(s) of . Only this class of distribution is unstable to strictly 
flute mode (ky = 0} propagation (Baldwin et al. , 1969). Pioneering work with 
this type of distribution was done by Tataronis and Crawford (1970) who investi- 
gated the stability of the distribution 

in an attempt to explain odd half harmonic emissions observed in laboratory 
devices. 

Shortly after the Ogo-5 magnetospheric observations, Fredericks (1971) 
seized upon the distribution Eq (2) (with a thermal spread in the perpendicular 
distribution) and was able to show the existence of unstable ”3/2” waves propa- 
gating within 10° of perpendicularity to the field for large but realistic values 
of cOp The electrons involved were thought to be ”hot” plasma sheet elec- 
trons: there was indication that ”3/2" emissions are most prevalent during 

substorm periods when the flux of plasma-sheet -like electrons at Ogo-5 alti- 
tudes is enhanced. 

Attempts to extent this work to more realistic (loss cone) distribution 
functions met in failure. Thus loss cone distributions of the Dory-Guest-Harris 
(DGH) (Dory et al . , 1965) type 
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(3) 
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were tried and shown to be electrostatically stable unless N became unrealis- 
tically large (>5) (Dory et al. , 1965). It was Young and his coworkers who 
first recognized the difficulty with small N; in those regions of w, k space where 
the imaginary component of Eq. (1) could be satisfied with 7=Imw>0 (corre- 
sponding to instability) the real component of Eq. (1) could not. However, the 
addition of a small amount of cold plasma to Eq, (3) alleviated the difficulty. 

The cold component was and is thought to be of ionospheric origin. The admix- 
ture which probably exists around LS 8 is optimum for Instability. 

The most recent work in the field, that of Ashour-Abdalla and Kennel 
(1976, cf, also references therein), also allows partial ”fill in” of the loss cone 
in order to keep growth rates acceptably small. (If growth is too fast, saturation 
occurs at a low emission amplitude, ) Ashour-Abdalla and Kennel thus analyze 
the stability of the distribution 


f(Ui, U|j) = 


1 




o;|| 

— II f 


.“iiH ^ “4 > “4 ' 


^ “l„ “IL ^ 


exp 


(4) 


This distribution is then the sum of three components: a cold bi-Maxwellian of 
density n^, a hot bi-Maxwellian of density, An^j, and a hot N = 1 DGH loss-cone 
distribution of density (1 - A) n^j. The ratio A/(l - A) between the densities 
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of hot Maxwellian and loss -cone plasma is referred to as the "fill-in" factor. 

The ratio jS is n^./n^ so tliat in Eq. (1) co^ is calculated using the hot plasma 
density n^ is then properly The Ashour-Abdalla and Kennel paper 

reports on a detailed parameter search for conditions which can produce "3/2" 
emissions and an analysis of the convective properties of the instabilities. It also 
contains a somewhat more detailed history of the evolution of theoretical ideas. 

Our work extends the plasma parameters search to the frequency regime 
f > 2fg. Since we are concerned with characterizing all categories of odd-half- 
harmonic emission discussed in Section II, we shall also he considering "3/2" 
emissions. Our results in this regime generally support the findings of Ashour- 
Abdalla and Kennel. 

We differ from previous work in one major respect: our hot electron dis- 
tribution parallel to B is a Lorentzian of order p; thus the distribution function 
which we analyze is that given by Eq. (4) with the substitution 

1 K2P-2)!! 

W»oi||^ " "h j-(2p-3)!!(V|f +<i|f^)P 

(0)!! = (-!)!! = 1 

The Lorentzian for several values of p is plotted in Figure 1 along with the 
Maxwellian for comparison. Note that as p becomes very large, the Lorentzian 
devolves to a 5-function at U|j = 0. We have chosen the Lorentzian for two 
reasons: 1) magnetospheric particle spectra frequently show a "power law" 
high energy tail (e.g. , the spectra shown in Anderson and Maeda, 1977) and 
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2) the ensuing complexity of the hot conti’ibution to Eq* (1) (cf Appendices) is 


greatly reduced when a Lorentzian is used. 

Let us write Eq. (1) symbolically in the form 


0 = 1 - 






( 5 ) 


where 


— ; andX|i*’= = 


D. 




We adopt the convention that k|[ >0, realizing that symmetric waves having 
kj[ <0 are equally possible. For the instabilities which we investigate in 
this paper Ajl > 2.5 and may range as high as 20. In Appendix A we show that 

H 


for large X, 



Here r(>/s) = /x + i'\/s Xyjj and erf(z) is the (complex) error function 


2 ^ 

erf(z) = — p=^/dt exp - 
0 


f2 


In addition, the validity of our asymptotic procedure in Appendix A requires 
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<9'(X||h ) ^ » these conditions being amply satisfied in all our 

computations, (Note that we can treat ii ~ X||jj ~ Xj^jj provided all are comparably 
large. ) 

The computational value of using a Lorentzian parallel distribution should 
be manifest in Eq. (6): neither sums nor integrations remain; the computation- 
ally most difficult component is the error function which is quickly and effi- 
ciently evaluated from a series representation; computer complex algebra 
readily separates Hjj into its real and imaginary parts and • 

In our investigations we have considered the p = 1 and p = 2 Lorentzians, 
for which 


H P=l = ■ 




l^J 


H 


H 


(7) 


where is the quantity in-f ->-in Eq. (6) (evaluated at s = 1 so that T-ii + iX||jj) 
and 


Were it necessary or desirable to study situations with yet higher p, one 
could continue this procedure and express completely in terms 

and its jit-derivatives. These derivatives may in tum be evaluated either 
numerically or algebraically. 
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The corresponding cold contribution 
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is of a form derived in several places in the literature (e. g. , Ashour-Abdalla 
et al. , 1975). In Eq. (9) is a modified Bessel Function of order n and Z 

is the well-lmown (Fried and Conte, 1961) plasma dispersion function 

+00 



exp - 

t- X 


While a strictly cold approximation to Eq. (9) is appropriate for some of our 
computations, X||^ and can in other cases be of ^(1) so that thermal correc- 
tions are important and several terms in the sum need be retained. Even under 
these conditions, however, the number of computations is fewer than for a ”hof' 
plasma. 

In solving Eq. (5) we consider only the situation close to marginal stability 
so that |)u.| <SC Denoting Hjj by Hv/e thus solve first the equation 

^2 

1 ( 10 ) 


in a manner outlined in Appendix B. Here is, in the usual sense, the real 

part of H, given ix to be itself real. We then determine a correction Aju:by 
solving the quadratic 
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( 11 ) 


^ H" + i = 0 

where in Eq. (11) the prime notation is nsed for differentiation and jUj. is that 
value determined from Eq. (10). Solutions to (11) are complex and thus pro- 
vide both fJL- and a correction to ii^ . The quadratic term in Eq. (11) is impor- 
tant only when = 0; when = 1, Eq. (11) reduces to the common- 

ly used prescription for calculating weak growth or damping. Having thus 
solved Eqs. (10) and (11) we check d posteriori that terms in (11) are in fact 
smaller than those in (10) and that other neglected terms are yet smaller. If 
this test of smallness fails, the roots are rejected. 
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IV. SOLUTIONS OF THE PLASMA DISPERSION RELATION 
Equation (5) is characterized by the following 6 parameters which can be 


OJ 


specified independently: , /3, A, T||^/T^jj , '^||c/T||h ■ would be an 

iTZ 

enormous task to examine a signtEicant volume of this 6-dimensional parameter 
space for instabilities, so we delimit the task somewhat by considering only the 
"isotropic" case where T|jjj “ Tj^h “"^h gyrofrequency 

scales all frequencies in the problem, and the gyroradius calculated using Q, and 
T^ is the fundamental unit of length. Our goal in this section is to summarize 
the gross effects of changing the remaining four parameters on instabilities in 
the frequency ranges described by our five proposed emission categories. In 
further discussion we consider the variables in the dispersion equation to be 
\h (calling it simply X) , and 6 = tan~^ (k^/kj|). 

Fig. 4 displays the frequency ju = cj/^2 of the fastest growing mode in each 
gyrohaz'monic band plotted against j3 = • The warm plasma density is constant 

= 20), T^/Tjj = 0.04, A = 0, and d = 86°. As n^, increases, (jl increases 
within each band, and higher frequency bands may appear. The dotted line indi- 
cates the cold upper hybrid frequency defined by 


— =^^uH = — 1(7:1“) 


( 12 ) 


The fact that instabilities do not occur for /it's much above is due to having 
Re’H dominated by the cold electrons for these values of Away from exact 
harmonics of ^2 , where is appropriate for instability, propagation in 
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the limit T^, 0 is possible only at = Muh* upper hybrid 

frequency has been investigated by Ashour-Abdalla and Kennel (1976) and Gaffey 
and LaQuey (1976). 

Figure 5 demonstrates the complicated effects of the cold plasma temper- 
ature on the range of the instability when > 2. Here |8 = 0.5, A = 0, and = 

3.32 (or = 20). The range 2.5 < A < 12 and 82" < $< 89“ has been sampled, 

and the region of corresponding unstable ju's is bounded by the solid curves shown. 
An asterisk indicates the fastest growing mode within each band at the corre- 
sponding T^/Tg value. For T^/T^ < 0. 01, note that (a) the lowest band is very 
narrow and lies close to [i~ 2 and (b) p is more nearly an absolute upper 
bound on fi than it is at higher /T^j . For T^, /Tj^ >0.01, bands above 4H and 
5S2 become unstable, and the instability width tends to narrow in the second and 
third bauds. The appearance of these higher bands is due to the increase in Ac 
with increasing T^. As A^ approaches''ll), thermal effects become important 
even for the "cold" electrons and several ln(A^/2) contribute to the sum in 
Hj, The result is that solutions with M>J“uH possible. 

In addition to shifting the unstable frequency range, the principal effect 
of increasing Tp/Tjj is to reduce most instability growth rates. Figure 6 plots 
the maximum growth rate as a function of T^/Tj^ for the parameters of 
Fig. 5, For T^^/Tj^ > 0. 04, further increases in T^ tend to stabilize the plasma 
in all bands. The lack of strong instabilities in the first, fourth, and fifth 
bands as T^,/Tjj 0 occurs because when T^ is very small, only modes near 
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can have large growth rates, The stabilizing effect of the cold plasma 
temperature arises primarily from the fact that <0 always. This 

gives a Landau damping effect which opposes theHjj contribution 

which drives the instability. 

For |3> 0. 05, the real part of the dispersion relation is dominated by 
the cold plasma component. Increasing Ujj should therefore increase the 
growth rate 7 without changing ix by much. This is verified by Fig, 7, which 
plots the most unstable frequency vs. j3 = with T^/Tj^ = 0. 10, A = 0, 

and - 5 = 2.45). The "error bars" delimit the frequency range in 

which 7 > Tjnax/^ hence give an estimate of the bandwidth. The solid, dashed 
and dotted lines correspond to T„,ax> 0-01 > 0-001 ^Tmax ^ 

respectively. Decreasing iijj tends to reduce as expected. The small vari- 
ation in [X which does exist is apparently an increase with decreasing n^^. 

We have already discussed the addition of a hot "fill-in" component which 
reduces the sharpness of the loss cone. Figure 8 shows the maximum growth 
rate 7^^^^^^^ in each band vs. A for = 3.32, n^, /n^ = 0. 5, and T^,/Tjj = 0. 04. 
Increasing the fill-in factor A decreases 7^^^^^^ in all bands but tends to affect the 
higher bands more strongly. In fact, for A - 0. 05, we found no unstable solu- 
tions above m = 5. This preferential stabilization of higher frequency bands 
when the loss cone is filled in is more pronounced for our Lorentzian parallel 
distribution than it would be for a Max^vellian of comparable temperature. The 
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fill-in component reduces growth rates because its contribution to i® 

always negative. 

Figures 4-8 reveal some of the complexities of the dispersion relation. 
The gross effects of changing each of the major plasma parameters are sum- 
marized below based on these Figures and an investigation of a much wider 
range of parameters. 

Cold plasma density . The cold plasma density n^ primarily determines 
the cold upper hybrid frequency and hence governs the number of unstable 
frequency bands. The highest unstable frequency always lies below ju-uh 
when Mjnax slightly above ja “jj if >2. If a particular band 

is unstable at some value of , the frequencies within that band will increase 
if n^ is raised further. 

V 

Hot plasma density . The hot plasma provides the free energy to drive 
the instability. Growth rates tend to increase with n^j , but the unstable fre- 
quency range is largely independent of n^j. 

Cold plasma temperature . The role of T^,/Tj^ is a complex one. As 
T^/T^ 0, large growth rates occur in the one or two bands closest to 

while other bands have much smaller growth rates. Low frequency bands 
(whose frequencies fx satisfy - [p]) > 1) occur near the top of the band and 
have small bandwidths. For T^/Tjj > 0. 01, one or two bands above may 
become unstable. At high cold plasma temperatures (0. 05 < T;./Tjj < 0. 2 ), the 
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unstable frequency range in each, band is shifted downward, and all growth rates 
are reduced. All instabilities disappear for Tj./T^ > 0. 2. 

Loss cone fill-in . Reducing the slope of (Uj^) reduces growth rates in all 
frequency bands but tends to suppress higher frequency bands much more 
severely than lower ones. 
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V. PLASMA MODELS FOR EACH EMISSION CLASS 


In the previous section, we examined the effects of changing various 
plasma parameters in our model on the solutions to the dispersion relation. We 
now proceed to use this information to infer plasma conditions which produce 
instabilities in the appropriate frequency ranges for each category of emission 
defined in Section II. The feet that we are successful in achieving this goal by 
variations in a single plasma model strongly supports the validity of the model. 
This single plasma model h5?pothesis is also suggested by the fact that there 
often seems to be a smooth transition in the Imp 6 data from one emission cate- 
gory to another {Shaw and Gurnett, 1975). We feel that the model is sufficiently 
successful that it can be used in a predictive fasliion to infer the unlmown cold 
plasma density and to a lesser extent the cold plasma temperature and loss cone 
fill-in factor whenever gyroharmonic emissions are observed. 

Using Section IV as a guide, we pick for each class of emissions a com- 
H 

bination of parameters — — , j ^ which produces frequencies and 

bandwidths characteristic of that category. By varying X and 6 , we determine 
the largest 7 corresponding to each unstable p and thus construct curves of 
maximum growth rate as a function of frequency. These curves should repro- 
duce the features of the emission class while having realistically small growth 
rates. Using Ashour-Abdalla and Kennel (1976) as a guide, we require temporal 
growth rates to satisfy 7^ 0. 01. Although we have not systematically calculated 
convective growth rates as was done by Ashour-Abdalla and Kennel, such low 
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temporal growth rates usually produce realistically small convective growth 
rates (k^aJSl < 1 where kj is the spatial growth rate). 

From Figures 4 and 5 and from Ashour-Ahdalla and Kennel (1976), the 
first two classes of emissions (”low 3/2” and ”3/2”) definitely require 
J^max ^ ^uH' somewhat less certainty we also infer that < 2 and 

0.05 < T /T„ < 0.2 when these Class 1 or Class 2 emissions are observed. 
Figure 9 represents 7 vs. ja for = 1. 34 (Curve 1) and = 1. 73 (Curve 2) 
with Tj./Tjj = 0. 10 and A = 0. Both curves show substantial growth rates above 
ju = 2 which can be suppressed completely by filling in the loss cone. The small 
curves (3 and 4) within the growth rate profiles represent the same parameter 
set with the fill in factor A = 0. 1. Hence, curve 3 represents the growth rate 
profile for a typical Class 1 (Ac^iax — emission while curve 4 represents a 
Class 2 (1.2 < Mmax — emission. Both curves have small growth rates and 
narrow bandwidths. The Class 2 or ”3/2" emission could be shifted closer to 
IX = 1. 5 by decreasing /T^^ or increasing . 

Figure 10 represents an attempt to replicate the features of "multiple 
half harmonic” or Class 3 emissions. The cold upper hybrid frequency [x^^ ~ 
3.32, A = 0, njn^ = 0,50, and T,,/Tjj =0,10. Higher cold plasma densities are re- 
quired than for ”3/2” emissions, and large fill in, which preferentially sup- 
presses the higher frequency modes, cannot be tolerated. The frequencies and 
bandwidths agree well with those shown in Fig. 1, and the successively smaller 
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growth rates as one goes to higher bands are in accord with the apparent mono- 
tonic decrease in intensity of the bands in Fig. 1. The unstable modes which 
occur just above ja = 5 propagate at a more nearly perpendicular angle to B 
(0 > 89®) while lower frequency modes tend to propagate most strongly in the 
range 80® < 0 < 86°. Modes with 6 ->• 90° have small Xj|jj and have frequencies 
very near the gyroharmonics. Instabilities in such modes tend to be stronger 
with our Lorentzian parallel distribution than with a Ma:xwellian. 

Since Class 4 o : diffuse band emissions almost always occur above 10 kHz, 
and since we have not yet examined any data above 10 IcHz, we have not presented 
a detailed parameterization of such emissions. However, preli m inary calcula- 
tions suggest that diffuse band emissions can also be explained by this same 
basic electron model, and we are continuing to investigate this possibility. 

Class 5 or f ~ f p narrow band gyroharmonic emissions can be achieved by 
this model only if >> n^j. This condition arises because the model cannot pro- 
duce instabilities with frequencies much above , and for cOp^ = £2^ , 

^UH <y CO + 1 Wpt/£2. In fact, if the cold upper hybrid frequency is 

exactly equal to the total plasma frequency - co ^/£2), one can easily show 

COpt 

that )3 = 1. It is which is determined from the lower cutoff of the 

trapped continuum radiation. If in Fig. 2, then the model requires 

Up/uji 80 since fp 9fg. However, we have seen that if > 0, frequencies 
somewhat above the cold upper hybrid can be generated and hence a smaller 
n^, /n^j ratio would be permitted. In some cases , gyroharmonic bands similar 
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to those in Fig. 2 are displaced just above fp; this situation requires a larger 
Uj, /Ufj than the one with gyroharmonic emissions at or just below f p. 

Figure 11 shows that gyroharmonic emissions at or above f p can be pro- 
duced if n^ /opj is sufficiently large. Since there is less free energy available 
from hot electrons to heat the cold plasma and since large cold plasma densities 
tend to reduce growth rates, the ratio T^,/Tjj is probably smaller for f ~ 
emissions than for ”3/2" emissions. In Fig. 11a, in which = 10, = 

0. 01, A = 0 and ~sj2\ ~ 4.58, a single narrow band emission near 
fp ~y/22 fg = 4.69 fg is the dominant feature. In Fig. lib, the cold plasma tem- 
perature has been raised to /T^ = 0. 02, An emission band appears above 
ju = 5, the band just below n = 4: might also be observable, while the band nearest 
fp has a lower growth rate than in Fig. 11a. This Figure suggests that (a) 
emissions which are clearly above and separated from fp, such as those in Shaw 
and Gurnett (1975, Fig. 4) are possible with this model and (b) the appearance of 
two or three bands near fp requires a higher T^/Ty than a single emission band. 
Multiple bands just below fp, such as those shown in Fig. 2, require a somewhat 
smaller nj,/njj than the value used in Fig. 11 (e.g. , 3 instead of 10). 

The results of this section are summarised in Table 2, which gives a 
reasonable estimate of the plasma conditions which we believe produce each 
type of emission. The ”3/2” model is in essential agreement with Ashour- 
Abdalla and Kennel (1976). This Table can be used as a crude diagnostic tool 
whenever gyroharmonic emissions are observed on a particular spacecraft. 
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TABLE 2 


Proposed Plasma Conditions for Each Class of Electrostatic Emissions 


Class 

Frequency Range 

=yi + 

Comments 

1. Low “3/2” 

1 ^ f/fg ^ 1.2 

7“ 1-5 

n^/uj^ « 1 (typically 0.1) 
0.05^T^/Tjj <0.15 

A < 0.1 (some fill in of loss cone 
distribution) 

2. “3/2” 

1.2gf/fg<2 


Essentially same as Class 1 but with 
larger cold plasma density 

3, “n+ 1/2” multiple 
half harmonics 

Wfg>2 


n,/nH <1,0.05 <T,/Th <0.15 
0 

4. “Diffuse band” 



n^/ng > 1 (?) 

T^/T.j -0.01 (?) 

Requires further study 

5. “f ~ fp” narrowband 


^uH _ 

n^/nH > U» 1 (?)) 
TJT^ <0.02, A «0 











VI. THE DIAGNOSTIC USE OF OUR PLASMA MODEL 


We have shown in the last section that each of our 5 categories of odd- 
half-Iiarmonic waves can be explained by suitable variations of parameters in a 
single plasma model. Table 2 summarizes the specific values or narrow ranges 
of n^/njj , T(,/T[j, and A characteristic of each class of emission. We suggest 
here that gyroharmonic wave observations be used in conjunction with Table 2 
and magnetometer measurements of the steady B-field as a diagnostic tool for 
determining local plasma parameters. We respond to the potential argument 
that the observed waves are not generated locally by pointing out that electro- 
static oscillations, unlike electromagnetic radiation, are rapidly damped as 
they propagate into regions which are not at least marginally stable to their 
growth. 

The plasma parameter which, is most definitively determined by observing 
odd-half ^harmonic emissions is /luy . In this section we shall develop further 
its relationship to emission frequencies, first by presenting in more detail re- 
sults of our theory and second by applying these theoretical results to observa- 
tions, both individual and statistical, made by the Imp-6 experiments. These 
observations are all consistent with our theory. 

While it is sometimes possible to deduce n^, by measuring n^ from the 
continuum cut-off and subtracting the warm plasma density as measured by 
particle detection experiments (Gurnett and Franit, 1974), this method ibecomes 
inaccurate when n^. n^j. The results of our calculations, previous 
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theoretical work (Young et al. , 1973; A shour-Ahdalla and Kennel, 1976), and 
the one observation presented by Gux'nett and Frank (1974) all suggest that 
n^ /njj 1 is a necessary condition for "3/2” emission to occur. Further- 
more, the appearance of these "3/2" waves is very common in the outer mag- 
netosphere. We conclude therefore that the use of the subtraction technique to 
calculate n^ is frequently a suspect one. 

Figure 12 represents our best extimate, based on solutions of the linear 
dispersion relation, fi^f the cold upper hybrid frequency as a function of the maxi- 
mum observed frequency of magneto spheric odd-half-harmonic emissions. The 
"error bars" are uncertainty estimates made under the assumption that the fill 
in factor A is <. 05 (which in turn assumes that cold plasma heating is the domi- 
nant stabilization mechanism (Ashour-Abdalla and Kennel, 1976; Grume, et al. , 
1972)). For Mmax s^lways so that we can deduce a very accurate 

lower limit for the cold plasma density. If in fact A >$> . 05 then the upper limit 
is higher, since under such circumstances may be -large and strong fill in 
can quench all but the "3/2" emissions. The dip at —1 reflects our 
computationally observed fact that for Classes 1 and 2, (^-iways), 

while this is not necessarily true if >2(Classes 3-5). We also observe 
computationally that A must be near zero for ^ ^ hence the error 
estimates in Fig. 12 are quite accimate for large . 

The estimates of { jl °^ are converted to densities in Fig. 13, using 
to 

— experimenters to estimate the cold 
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plasma density whenever gyroharmonic bands are observed. When the con- 
tinuum cutoff is simultaneously observed, one can estimate n^ = n^ ~ n^,. It is 
usually possible to estimate to within 5% using the electric field data and 
the measured local magnetic field strength. 

We have previously stated that Class 5 (f ~ f p gyroharmonic) emissions 
with fp fg probably require n^. » Uj^. Whereas Classes 1-3 emission tend 
to occur during geomagnetically disturbed times when n /n„ is expected to be 
small in the outer magnetosphere (cf. Anderson & Maeda, 1977), Classes 4 
and 5 are most common during quiet times (Shaw and Gurnett, 197 5) when this 
ratio is more likely to be large , as predicted by our theory. We predict that when- 
ever f fp gyroharmonics are observed, n^, is anomalously Ingh. At such times 
one would expect that particle detection experiments which measure particles only 
above some minimmn energy threshold can give a significantly lower electron 
density than tlie density measured from the trapped continuum lower cut-off. We 
are currently investigating the data for verification of this hypothesis. 

Our discussion of the variation of nj./nj^ from class to class is pretty much 
independent of the precise form of the electron distribution fimction, so long as 
it has a cold component and is loss cone unstable. Such is not the case with T /T„ 
and Table 2 should be used vifch according caution. The hot distribution fimction, 
which we have used for our calculations has a loss cone which is unrealistically 
large for the real outer magnetosphere, A more plausible distribution with a 
peak in fj^j^ at a much smaller value of Uj^ (relative to ) can likewise be 
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unstable (Young et al. , 1973). The characteristics of the instability depend 

strongly on the positive slope (9fjL/avj^ > 0) portion of the distribution function 

and are virtually independent of the tail (9fx/3vjL ^ 0) so long as it is "gentle. " 

In appl 3 dng Table 2, T^j should be interpreted as the effective temperature 
mv ^ 

Tr eff corresponding to the speed at which fj^ maximizes, since the 

important low region can generally be fit by a DGH distribution of the type 
we use. An important consequence is that given a hot electron component with 
a 1 keV thermal spread typical of the plasma sheet, Tr can be 100 eV or 
smaller, and hence from Table 2 T^. 10 eV, a realistic value for plasma of 

plasmaspheric or ionospheric origin. 


The use of Fig, 13 and Table 2 as a diagnostic tool can be illustrated using 


Figures 1 and 2. The emissions at 0715 of Fig. 1 are Class 3 (multiple half 
harmonic emissions) with = 6.1 and = 1. 15 kHz. (All B values are 
obtained from the GSFC magnetometer data,) Hence, we estimate 
co^jQ.^ «s30 fi-om Fig. 13, and therefore n^ ^ 0.5 cm"^ . From Table 2, A 0 
and 0. 05 ^ T^,/Tr ^ 0.2. Trapped continuum radiation is not observed, so the 
total electron density is not available. Figure 2 is a Class 5 emission with 
^max " ^ 0537. Using Fig. 12 as 

a guide, we estimate ju“r ^ S, so n^ 0.38 cm“^. Since the total plasma density 
(from fp) is 0. 50 cm"^ , the ratio n^,/nR 0.38/0.12 ^3.2. When the den^ 

sity ratio predicted by this method is quite sensitive to the estimate of 
From Table 2, we estimate TjT^ ~ 0.02 and A 0. 
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Figui’es 12 and 13 are also supported by a statistical survey of the data 
taken from a number of orbits. Figure 14 shows as a function of f^ the highest 
value of f/f observed during 15 minute intervals from orbits 178 through 183 

a 

of Imp-6 whenever odd-half -harmonic waves are present. The data were all 
taken on the day side of the outer magnetosphere at moderate magnetic latitudes 
(20° ^ iXjjj I < 55°). The tendency to see higher relative frequencies when the g 3 o;o- 
frequency is low can be understood if the cold plasma density changes much more 
slowly than f^ once the spacecraft is more than ~lE,j 7 beyond the nlasmapause. 

In that case, and hence will tend to increase as f^ decreases, and the 

theoi'y predicts that higher values of y (and hence more bands) may be observed. 

The data used in Fig. 14 also suggests thatAtjnax decreases at high magnetic 
latitudes . We divided the data into high (IX^^j | ^40°) and low (iX^j^ I <40°) mag- 
netic latitude cases. 82% (23 of 28) of the high latitude observations were Class 
1 (Mn,ax emissions. However, only 53% (29 of 55} of the low latitude ob- 

servations were Class 1, and all of the Class 3 (m > 2) observations were at 
low latitudes. This effect can be theoretically interpreted using our model. At 
constant radial distance R/Rg , f^ increases and n^, is probably constant or de- 
creasing at high latitudes. Hence and thus would tend to decrease. 

This probably accounts for the preponderance of Class 1 < 1-2) observa- 

tions in the Imp -6 data. Explorer 45, which always remains at low latitudes, 
observed more of the traditional y « 3/2 emissions (R. R. Anderson, private 
G ommunication) . 
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VII. CONCLUSIONS 


We have conducted a theoretical analysis of the banded electrostatic 
emissions which occur between consecutive harmonics of f in the outer mag- 
netosphere and have compared our theoretical predictions with Imp 6 electric 
and magnetic field data. Our results are summarized below. 

First, we have constructed a scheme for classifying these emissions 
(Table 1) . Our five categories include the widely studied odd-half -harmonic 
emissions as well as the waves reported by Shaw and Gnrnett (1975). Second, 
we have extended successfully the linear plasma model of Young, et al. (1973) 
and Ashour-Abdalla and Kennel (1976), which can explain "3/2" emissions, to 
the other categories of emissions. Our results indicate that a single plasma 
model consisting of a hot loss cone electron distribution and a cold electron 
component can generate instabilities in the frequency ranges appropriate to 
each category. Hence, all of these emissions appear to arise from the same 
physical mechanism. Figure 13 and Table 2 can be used as diagnostic guides 
for estimating the cold plasma density and temperature whenever banded electro- 
static emissions are detected. Finally, our results indicate that certain types 
of emissions require the cold plasma density n^ to exceed the hot plasma density 
n^j. Hence, although the hot loss cone component usually dominates the density 
in the outer magnetosphere, there are times when n^, njj. 
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0"I0 kHz 


B-G77-6I 


lMP-6 

U OF IOWA PLASMA WAVE EXPERIMENT 

ORBIT 178 MARCH 19, 1973 



fg = 1.15 kHz 


Figure 1. An example of multiple half harmonic emissions observed by IMP 6 at S « DKeou the day 
side of the magnetosphere. As many as seven g 3 n:oharmonic bands are simultaneously 
observed. The intense ’'3/2” emission which occurs between 1 and 2 kHz at 0700 is 
almost continuously observed from 0500 to 0800. The emissions below 1 IcHz are electro- 
magnetic and are not considered in this paper. 
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FREQUENCY 
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U OF IOWA PLASMA WAVE EXPERIMENT 

ORBIT 84 FEBRUARY 21, 1972 
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fg = 0,70 kHz 
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Figure 2. An example of gyroharmonic emissions occurring near the electron plasma frequency fp . 

The shaded region is trapped continuum radiation whose lower cutoff is fp . Up to three 
gyroharmonic bands are observed. 
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Figure 3. A comparison between a Maxwellian parallel velocity 

distribution and the corresponding p**’ order Lorentzian 
distribution. The p = 2 Lorentzian has a similar shape 
for V|| ^ a|| but has more high energy particles than 
does tile Maxwellian. 
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Figure 4. The dependence of the frequency ^ = f/fg of the fastest growing mode 
within each gyroharmonic band on the cold plasma density. The 

parameters ~ 20, T^/Th = 0.04, 6 - 86°, and A = 0. As 

Uc/nn increases m increases within each band. The highest frequency 
tracks the cold upper hybrid frequency {dotted line). 



Figure 5. The dependence of range of the instability within each band on the 
cold plasma temperature for j3 = 0. 5, A = 0, and mSh = 3.32. The 
region of unstable frequencies within each band is enclosed by the 
solid curves, and an asterisk indicates the fastest growing mode 
within each band as a function of T^j/T^. The temperature ratio has 
a considerable effect on the range of the instability. 



Pigure 6. The maximum growth rate within each of the five gyroharmonic 
bands shown in Fig* 5, plotted against the temperature ratio Tq/Th. 
For Tc/Th < 0, 01, strong growth is seen only in the second and 
third bands which are closest to = 3* 32. For Tc/Th > 0. 04, 

Tmax decreases monotonicaily with Tc in all bands. 
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Figure 7. The effect of the hot plasma density on the range of unstable fre- 
quencies. The most unstable frequency is plotted vs. Uj./nj^ with 
Tc/Th = 0. 10, A = 0, and Uj. fixed ’'error bars" 

estimate the bandwidth based on the frequency range for wMeh the 
maximum growth rate 7max < The solid, dashed, and dotted 

curves correspond to 7max > 0 . 01 , 0. 001 < Tniax < 0 . 01 , and 

7max 0. 001 respectively, hicreasing tends to increase 
while leaving frequencies relatively unaffected. 
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Figure 8. A plot of vs. fill in factor A in each of tlie five bands wMch 

occur when Tc/Th = 0. 04, = 3. 32, and n^/n^ = 0.5. Growth 

rates in all bands decrease monotonically with A , but higher fre- 
quency bands are suppressed most severely. 
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Growth rate 7 vs. frequency fi for ;^,,h = 1. 34, |3 =0.1 (Curve 1) and 
Muh = 1.73, 13 = 0.2 (Curve 2) with T^./Th = 0. 10 and A == 0. Unstable 
waves above M = 2 can be suppressed by setting loss cone fill in 
factor A = 0 . 10 (Curves 3 and 4). Hence Curves 3 and 4 reproduce 
the features of Class 1 (low 3/2) and Class 2 (3/2) emissions. 
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Figure 10. 7 vs. ii for = 3. 32, A= 0, p = 0, 5, and Tp/Tjj = 0. 10, The 

frequencies and bandwidths are characteristic of Class 3 (multiple 
half harmonic) emissions. 



Figure 11. An example of gyroliarmonic emissions near the plasma frequency. 

Fig. 11a gives 7 vs. /u for |3 = 10, A = 0 , Muh = 4.58, fp/fg = 4.69, 
and T(,/Ti{ = 0. 01, The dominant feature is a single strong emission 
near f p. In Fig. 11b, T(./Th = 0. 02, and an emission band above 
5fg would be seen. Class 5 emissions include both single band and 
multiple band emissions near fn. 





Figure 12. An estimate, based on solutions to the linear dispersion relation, of 
the cold upper hybrid frequency jLigjj vs. maximum observed fre- 
quency Mniax • "error bars" are. uncertainty estimates assuming 
A <0.05. This figure can be used to estimate mSh from electric and 
magnetic field data whenever emissions of any class are observed. 




Figure 13. The estimates of Fig, 12 converted to (or equivalently n^,) 

vs, maximum observed frequency • This figure can be used as 
a diagnostic tool to estimate the cold plasma density nc . 
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Pigure 14, A plot of 83 individual electric field observations talcen at 15 minute 
intervals from Imp 6 Orbits 178-183 whenever "odd-haLf-harmonic" 
emissions are seen. The graph shows the maximum observed 
frequency Mmax (hi units of fg) vs. fg. Higher relative frequencies 
tend to occur when the magnetic field is weakest, and almost all 
emissions between fg and 2fg occur below 1, 5^ . The region above 
the solid curve corresponds to frequencies above 10 IdHz and cannot 
be observed in the data. 
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APPENDIX A 


EVALUATION OF HOT ELECTRON CONTRIBUTION 
TO THE DISPERSION RELATION 


Since f = the hot and cold electron components contribute sepa- 

rately to the dispersion equation, Eq. (1). We consider here the hot contri- 
bution 
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where the -notation is that introduced in conjunction with Eq. (5) , and we 
have partitioned 
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where the Ij^'s are modified Bessel Functions of order n and the notation 
been used. In a similar manner 
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For the parallel integration we insert the Lorentzian form and define 
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The first U|| -integral in (1) thus becomes 
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where N denotes the normalization factor 
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the latter form of (A- 5) following by contour integration-closing in the negative 
imaginary x~plane and picking up the contribution from the simple pole at 
x = -i-v/s7 (According to our convention, integration in the complex co-plane is 
to be carried out with Im to > 0, so that for integration purposes > 0.) 

In exactly the same way we evaluate 
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If one of the differentiations in (A-5) is carried out explicitly, Hjj can be ex- 
pressed in the particularly convenient form 
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If we define further 
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In. Eq. (A-9) 
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and use has been made of the identity 
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We now proceed to evaluate the summation 
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in the limit that z -> <» . Note that 
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This result is most easily proven by working in the reverse direction, i. e. , 
expanding exp(z cos 0) in terms of a Bessel Series and performing then the 
0-integration. Further 
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The form (A-12) is ideal for asymptotic analysis; for large z, the major 
contribution to the integral comes from a small region around 0 = 0 where 
cos 0 = 1, and we may series expand the cosine function 
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Convergence in (A-13) comes from the exp - factor and has essentially 
occurred by the time 0 ^ term is then an O' correction to the 

dominant term in {A-13). The&izdfi) correction is of and still smaller. 

wJ 

Extension of the upper integration limit to introduces an error exponentially 
small in z. We retain the dominant and first non-dominant term in g and 
write 
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where the error function has its usual definition 
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Carrying out the differentiation indicated in (A-14), we obtain 
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Finally we conclude 
+00 

^ r-n \ 2z } V2z 


(A-15) 


sin 2?rr 


I - cos ttP 


- i erf 


Vv^yJL 24 zI^ z 


+ Sir 

24z2 


(=-f) 


(A-16) 


A-6 


The leading order term is the one of order The asymptotic smallness of 

corrections depends on P^/z remaining of ©"(I) or smaller as 
For P^/z""0'(l) all correction terms in (A-16) are of comparable magnitude 
and are small compared to the dominant by the anticipated l/z. 

We now incorporate the leading order term of Eq. (A-16) into Eq. (A -9). 
While it would appear that terms in (A-9) having as a multiplicative factor 


are asymptotically dominant, operation by 


H) 


drops their magnitude by 


one compensating order in . We are justified in retaining only the lowest 
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order term in (A-16) because neither the operations 7 —;= nor 

o^/s 


changes 

the relative ordering of the dominant and correction terms. Our final result is 
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APPENDIX B 


OUTLINE OF COMPUTER SOLUTION TO THE DISPERSION RELATION 


We here discuss computational detail briefly beginning with a summary of 
our treatment of We choose initially values of ’ ®iH ^®|lH > propa- 

gation direction 0 = tan ~^ Our work thus far has concentrated on the ”iso- 

Kii 


tropic" case discuss only these results. Given the values of 

these three chosen parameters , is determined. We next vary Mj., typically 
between = 1. 02 and ju^ = 5. 98, in steps of . 04, evaltmting and 

using the computer to sort real and imaginary parts. In the process we also 
compute the first three Mr derivatives of both the real and imaginary parts at 
each fx^ value. We store these values of and its derivatives on tape, tabu- 
lating sep.’rately the contributions of the "loss cone" and "filled in" components, 
which are characterized respectively by the multiplicative factors (1 - A) and 
A in Eq, (6). This table is then repeated for various pairs of Xj^jj and 6, Xj^jj 
typically varying between 2. 5 and 16 in a variable step size (.5 < AX ^ 2) and 9 
taking on the values 89, 5®, 88° , 86°, 82°, and 76°. 

We combine each tabulation of characterized by fixed Xj^j^, 0, and^r* 
with several "realizations" of cold plasma. In characterizing the cold plasma 
we still have the freedom to choose /? and ctja^ . (Again, we have as yet con- 
sidered only «[|c = aj_c)- Once these parameters are chosen H^^rcai ^c.imag 
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are calculated from Eq. (9). This procedure is a standard one, and further 
discussion of it will be omitted. 

Finally is evaluated from Eq. (10) and fx. and are deter- 

mined from Eq. (11). [Apj. is typically <0,02 for accepted results.] 

While the sequential events of computation are as detailed above, para- 
metric dependencies are ultimately determined by cross referencing. For 
example, one would determine the dependence of and /u. on Xjlh hence 
^ic ) for fixed S,j3, and by first generating several tables each character- 
ized by a different value of , Each table would have as "output" the values 
of CO p /i2^, /ij, and (Ap)j. as a function of the "input" One would then further 
search from table to table at fixed cOp^j/Sl^ to determine, usually by interpola- 
tion jjL^ and + (A}jl)^ as functions of . 
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